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^ ■ Abstract 

^ , The definition of Q-conditional symmetry for one PDE is correctly generalized to a 

' special case of systems of PDEs and involutive families of operators. The notion 

. of equivalence of Q-conditional symmetries under a group of local transformation is 

\ introduced. Using this notion, all possible single Q-conditional symmetry operators 

' are classified for the n-dimensional (n ^ 2) linear heat equation and for the Euler 

ly-^ ■ equations describing the motion of an incompressible ideal fluid. 

o : 

O ■ The concept of Q-conditional symmetry called also nonclassical symmetry was introduced 
^ . by Bluman and Cole in 1969. This year is the year of the 30*^ anniversary of appearance 
CN ! of their pioneering paper Although the concept of Q-conditional symmetry exists for a 

^ \ long time and has various applications many problem of its theory are not solved so far. 

Before 1986 the nonclassical symmetry was only mentioned in few papers, e.g., in 0. 
I ■ The intensive application of Q-conditional symmetries to finding exact solutions of partial 
■ differential equations (PDEs) and the parallel search for their foundations was begun after 
^ ! publication of the papers of Olver and Rosenau in 1986 and 1987 ^ as well as the paper 
\ of Fushchych and Tsyfra in 1987 0. 

I The first correct definition of a Q-conditional symmetry operator for one PDE was given 

^ ■ in lH- Later it was generalized to involutive families of operators [6-8]. We stress that it can 
. be directly extended only to some special cases of Q-conditional invariance for systems of 

PDEs. For all the other cases this definition must be essentially modified and is much more 

complicated. 

In this paper we correctly generalize the definition of Q-conditional symmetry [6-8] to a 
special case of systems of PDEs and involutive families of operators. Further, we introduce 
the notion of equivalence of Q-conditional symmetries under a group of local transforma- 
tion. Using this notion, we can, first, classify all the possible Q-conditional symmetries and, 
correspondingly, all the possible reductions of systems of PDEs |]^, § and, secondly, essen- 
tially simplify the procedure of finding Q-conditional symmetries in some cases when the Lie 
symmetry group is sufficiently wide. 

Consider a system of k PDEs of the order r for m unknown functions u = [u^, . . . 
depending on n independent variables x = (xi, . . . , Xn) of the form 

L{x,U(^r){x)) = 0, L = (L\...,L^). (1) 

Here the order of a system is the order of the major partial derivative appearing in the 
system. The symbol U(^r) denotes for the set of partial derivatives of the functions u of the 
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orders from to r. Within the local approach system (|I]) is treated as a system of algebraic 
equations in the jet space J'-''^ of the order r. 

Consider also an involutive family Q of / differential operators 

= C{x,u)d^^ + Tl'\x,u)du^, where / ^ n, i^nk\\C{x,u)\ \ = I. (2) 

The requirement of involution means for the family Q that the commutator of any pair of 
operators from Q belongs to the span of Q over the ring of smooth functions of the variables x 
and M, i.e. 

Vs,p ^r^' = C'"\x.n): [Q%Q^]=C'"'Q''. (3) 

Here and below the indices a and h run from 1 to m, the indices i and j run from 1 to n, the 
indices s and p run from 1 to /, and the indices /x and v run from 1 to n — /. The sumation 
is imposed over the repeated indices. Subscripts of functions denote differentiation with 
respect to the corresponding variables. 

If operators (0) form an involutive family, then the family Q of differential operators 

Qs^ypQP^ where X'^ = \'p{x,u), det | |A"p| | ^ 0, (4) 

is also involutive. And family (^) is called equivalent to family (|^) [6-8]. 
Notation: Q = {Q'} ~ Q = {Q'}. 

By the Frobenius theorem, condition (|^) is sufficient for the system of PDEs 

Q^K]:=r^^'^(x,«)-r(x,«)— = (5) 
to be compatible. 

Denote the manifold defined by the system of algebraic equations L = in J^"^^ by C and 
the manifold corresponding to the set of all the differential consequences of the system of 
PDEs (D in J(^) by M: 

C = {(X, 'U(r)) G J^^'^ I L{x,U{r)) = 0}, 

M= {{x, U(^r)) e J(^) I • • • ^""<5'M = 0, fti G N U {0}, |a|: = ai + ■ ■ ■ + a„ < r}, 
where Di = dx^ + X] ""a i^u% is the operator of total differentiation with respect to the 



variable Xj, and u'^ ^ denote the variables in J'^'"\ corresponding to derivatives ^^^^ g^a„ 

Let the system L|;vi = do not includes equations which are differential consequences of 
other its equations. Moreover, let all the differential consequences of the system L\j^ = 0, 
the orders of which (as equations) are less than or equal to its order, vanish on £ ft A^. 

Definition 1. System of smaller PDEs is called Q-conditional invariant with respect to 
involutive family of differential operators ^ if the relation 

(QUL) = (6) 



holds true. Here the symbol Q^-^ denotes the rth prolongation of the operator 



\a\^r 

Denote the set of involutive families of I operators of Q-conditional symmetry of system ([1|) 
as B{C, I): 



B{C,l) = lQ = {Q\...,Q'} 



the system L = is Q-conditionally 
invariant with respect to Q 

Lemma [6—8]. Let system of PDEs 1^ be Q- conditionally invariant with respect to involu- 
tive family of operators ([^. Then, it is Q- conditionally invariant with respect to an arbitrary 
family of the form i. e. 

QeB{CJ), Q^Q =^ Q e B{CJ). 

An important consequence of the lemma is that we can study Q-conditionally invariance 
up to equivalence relation (^) which is defined on the set of involutive families of I operators 
as well as in B{C, I). Then it is possible for an arbitrary family of operators (Q) to choose the 
functions X^p{x,u) and, if it is necessary, to change enumeration of the variables xi, . . . , Xn 
in such a way that operators (^) take the following form: = dx^ + ^'^'^^^ dx^^^ + fj^'^du'^ . 
Operators generate a commutative Lie algebra. 

Let A{C) and G{C) denote the maximal Lie invariance algebra of system (|1|) and its max- 
imal local symmetry correspondingly. Now we strengthen the equivalence relation in B{C, I), 
given by formula (^, by means of generalizing equivalence of /-dimensional subalgebras of 
the algebra A{C) under the adjoint representation of the group G{C) in A{C). 

We use the following lemma for this generalization. 

Lemma. Let g be an arbitrary local transformation from G{C). Then the adjoint action of 
g in the set of differential operators generate a one-to-one mapping from B{C, I) into itself. 

Let Q = {Q'^} and Q = {Q'"^} be involutive families of differential operators. 

Definition. The families Q and Q are called equivalent with respect to a group G of local 
transformations if there exists a local transformation g from G for which the families Q and 
Ad{g)Q are equivalent. 

Notation: Q Q mod G. 

Definition. The families Q and Q are called equivalent with respect to a Lie algebra A of 
differential operators if they are equivalent with respect to the one-parametric group generated 
by an operator from A. 

Notation: Q Q mod A. 

Therefore, 

Q-gmodG 44 3geG:Q^Ad{g)Q. (7) 

Qr^Q mod A 44 3VeA: Q r^Q mod {e'^, eeU{0, 6) CR} (8) 
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Lemma. Formulas (0j and define equivalence relations in the set of involutive families 
of I differential operators. Moreover, if G is a subgroup of G{C) and A is a subalgebra of 
A{C) then formulas Q) and (Qj define equivalence relations in B{C,l). 

Comsider two examples. 

Example 1. Investigate Q-conditional invariance of the linear n- dimensional heat equation 

Ut = Uaa, where u = u(t,x), t = xq, x = (xi, . . . , x„), (9) 

with respect to a single operator (Z = 1). 

It is just the problem with n = 1 for which Bluman and Cole introduced the concept of 
nonclassical symmetry. In the one- dimensional case the problem was completely solved in 
IP]. That is why we pay our attention to the multidimensional problem. 

Lie symmetry of equation (^) is well known. In the one- dimensional case it was investigat- 
ed by Lie. The maximal Lie invariance algebra A(LHE) of @ is generated by the following 
operators: 

dt = d/dt, da = d/dxa, D = 2tdt + Xada, Ga = tda - ^XaUdu, I = udu, ^^^^ 

Jab = Xadb - Xbda {a < b), U = At"^ + AtXada - {XaXa + 2t)udu, f{t,x)du, 

where / = f(t,x) is an arbitrary solution of (^. 

Theorem 1. For any operator Q of Q- conditional symmetry of equation one of three 
following conditions holds: 

1. Q^Q\ where QOgA(LHE); 

2. Q = dn + Qng'^udu mod ASO{n) + A°°(LHE), where g = g{t,Xn) (gn ^ 0) is a 
solution of the one- dimensional heat equation, that is, gt = gnni 

3. Q = Ji2 + (pi9)udu mod AG{l,n) + A°°(LIIE), where ip = ip{9) is a solution of 
the equation (fgo + 2(f(f0 = 0, (fg ^ 0, 6 is the polar angle in the plane OX1X2. 

Here 

A^(LHE) = {f{t,x)du\f = f{t,x):ft = faa), 

AG{l,n) = {dt,da,Ga,Jab), ASO{n) = {Jab). 

It follows from Theorem 1 that there exist only three classes of the possible reductions on 
one independent variable for the linear multidimensional heat equation. 
The first class is formed by Lie reductions. 

The second class involves reductions which are similar to separation of variables in the 
Cartesian coordinates: 

u = g(t, Xn)v{uo, Un-i), where ljq = t, uJi = xf, (|) =^ vq = Vu. 

The third class is formed by reductions which are similar to separation of variables in the 
cylindrical coordinates: 

u = exp( f ip{0)d6)v{(jjQ, . . . , ujn-i), where uq = t, Ui = r, ojs = Xs+i, s = 2, n — 1; 
(D =^ Vo = Vn + u^^vi - Xuj^^v + Vss- 
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Here A = —ipe — ip^ = const, (r, 6*) are the polar coordinates in the plane OX1X2. As the 
equation ipgo + 2LpLpo = has four essentially different (under translations with respect to 
6) families of solutions with (fg 7^ 0, there are four inequivalent cases for the third class of 
reductions (x 7^ 0): 

a) if = — K ta.n x6 : u = v{ujq, . . . ,ujn-i) cos x6, A = x^; 

b) <^ = xtanhx^: u = v{uJo, . . . ,ujn-i) cosh x6, A = — x^; 

c) y9 = xcothx^: u = v{u!o, . . . ,u!n-i) sinh >i6, A = — x^; 

d) = : u = v{ujq, . . . ,Un^i)6, A = 0. 

Example 2. Consider the Euler equations 

ut + iu-V)u + Vp = 6, divu = (11) 

describing the motion of an incompressible ideal fluid. In the following u = {u"'{t, x)} denotes 
the velocity of the fluid, p = p{t,x) denotes the pressure, n = 3, x = {xa}, dt = d/dt, 

—t — » — » 

da = d/dxa, V = {da}, A = V ■ V is the Laplacian. The fluid density is set equal to unity. 

Lie symmetry of system ([llD was investigated by Buchnev [|l^, |lT|. The maximal Lie 
invariance algebra A{E) of (|IT]) is infinite dimensional and generated by the following oper- 
ators: 

dt, Jab = Xadb - Xbda + u^d^b - u^'dua (a<6), 

= tdt - u^dua - 2pdp, D"" = Xada + u^dua + 2pdp, (12) 
R{m) = m%t)da + m^{t)d^a - m^^{t)xadp, Z{x) = x{t)dp, 

where = m°'{t) and x = xit) ^-^e arbitrary smooth functions of t (for example, from 
C°°{{tQ,ti),M.)). Let us investigate Q-conditional symmetry of ( |Tll) with respect to alone 
differential operator Q = ^°(t, x, u, p)dt + ^""{t, x, u, p)da + T]°-{t, x, u, p)du'^ + rp{t, x, u, p)dp. 



Theorem 2. Any operator Q of Q -conditional symmetry of the Euler equations (jll[) either 
is equivalent to a Lie symmetry operator of ^Tl\ ) or is equivalent (mod A{E) ) to the operator 

Q = d3 + C{t, X3, u^) du3 + x{t)x3dp, (13) 

where O ^ 0, Ca + CO = 0,Ct + (u^ + XX3) O + (C)' + X = 0. 

It follows from Theorem 2 that there exist two classes of the possible reductions w.r.t. 
independent variable for the Euler equations, namely, the Lie reductions and the reductions 
corresponding to operators of form (pIS)). 

Lie reductions of the Euler equations ([TT|) are investigated in [12-14]. 

An ansatz constructed with the operator Q has the following form: 

= y-^, = , V? = x^v^ + ip (t, f ^) , p = q + ^x{t)xl, 
where v"" = v°'{t, xi, X2), q = q{t, Xi, X2), ip = "^(^j '^^) is a solution of the equation 
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Substituting this ansatz into (|n|), we obtain the corresponding reduced system = 
1,2): 

vi + vH) + gi = 0, v^v] + [f) V X = 0, + t;^ = 0. 

The analogous problem for the Navier-Stokes equations 

Mi + (m ■ V)m + Vp- z/Am = 0, divM = [v^^) (14) 

describing the motion of an incompressible viscous fluid was solved by Ludlow, Clarkson, 
and Bassom in Their result can be reformulated as follows. 

Theorem 3. Any (real) operator Q of Q- conditional symmetry of the Navier-Stokes equa- 
tions is equivalent to a Lie symmetry operator of (|7^. 

Therefore, all the possible reductions of the Navier-Stokes equations w.r.t. independent 
variable are exhausted by the Lie reductions. Lie symmetry of system (|14D was studied by 
Danilov |]16|, [l^ . The maximal Lie invariance algebra of the Navier-Stokes equations (|14| is 
similar to one of the Euler equations (see (|T^)): 

A(NS) = {dt, Jat, D' + \D\ R{m{t)), Z{Cm- 

The Lie reductions of the Navier-Stokes equations were completely described in 
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